Local Dialogues Between the Endocrine and Exocrine Cells in the Pancreas.

For many years, it has been taught in medical textbooks that the endocrine and exocrine parts of the pancreas have separate blood supplies that do not mix. Therefore, they have been studied by different scientific communities, and patients with pancreatic disorders are treated by physicians in different medical disciplines, where endocrine and exocrine function are the focus of endocrinologists and gastroenterologists, respectively. The conventional model that every islet in each pancreatic lobule receives a dedicated arterial blood supply was first proposed in 1932, and it has been inherited to date. Recently, in vivo intravital recording of red blood cell flow in mouse islets as well as in situ structural analysis of 3D pancreatic vasculature from hundreds of islets provided evidence for preferentially integrated pancreatic blood flow in six mammalian species. The majority of islets have no association with the arteriole, and there is bidirectional blood exchange between the two segments. Such vascularization may allow an entire downstream region of islets and acinar cells to be simultaneously exposed to a topologically and temporally specific plasma content, which could underlie an adaptive sensory function as well as common pathogeneses of both portions of the organ in pancreatic diseases, including diabetes.

The effect of forskolin and the role of Epac2A during activation, activity, and deactivation of beta cell networks.

Beta cells couple stimulation by glucose with insulin secretion and impairments in this coupling play a central role in diabetes mellitus. Cyclic adenosine monophosphate (cAMP) amplifies stimulus-secretion coupling via protein kinase A and guanine nucleotide exchange protein 2 (Epac2A). With the present research, we aimed to clarify the influence of cAMP-elevating diterpene forskolin on cytoplasmic calcium dynamics and intercellular network activity, which are two of the crucial elements of normal beta cell stimulus-secretion coupling, and the role of Epac2A under normal and stimulated conditions. To this end, we performed functional multicellular calcium imaging of beta cells in mouse pancreas tissue slices after stimulation with glucose and forskolin in wild-type and Epac2A knock-out mice. Forskolin evoked calcium signals in otherwise substimulatory glucose and beta cells from Epac2A knock-out mice displayed a faster activation. During the plateau phase, beta cells from Epac2A knock-out mice displayed a slightly higher active time in response to glucose compared with wild-type littermates, and stimulation with forskolin increased the active time via an increase in oscillation frequency and a decrease in oscillation duration in both Epac2A knock-out and wild-type mice. Functional network properties during stimulation with glucose did not differ in Epac2A knock-out mice, but the presence of Epac2A was crucial for the protective effect of stimulation with forskolin in preventing a decline in beta cell functional connectivity with time. Finally, stimulation with forskolin prolonged beta cell activity during deactivation, especially in Epac2A knock-out mice.

Stromal Interaction Molecule 1 Maintains ß-Cell Identity and Function in Female Mice Through Preservation of G-Protein-Coupled Estrogen Receptor 1 Signaling.

Altered endoplasmic reticulum (ER) Ca2+ signaling has been linked with ß-cell dysfunction and diabetes development. Store-operated Ca2+ entry replenishes ER Ca2+ through reversible gating of plasma membrane Ca2+ channels by the ER Ca2+ sensor, stromal interaction molecule 1 (STIM1). For characterization of the in vivo impact of STIM1 loss, mice with ß-cell-specific STIM1 deletion (STIM1Δß mice) were generated and challenged with high-fat diet. Interestingly, ß-cell dysfunction was observed in female, but not male, mice. Female STIM1Δß mice displayed reductions in ß-cell mass, a concomitant increase in α-cell mass, and reduced expression of markers of ß-cell maturity, including MafA and UCN3. Consistent with these findings, STIM1 expression was inversely correlated with HbA1c levels in islets from female, but not male, human organ donors. Mechanistic assays demonstrated that the sexually dimorphic phenotype observed in STIM1Δß mice was due, in part, to loss of signaling through the noncanonical 17-ß estradiol receptor (GPER1), as GPER1 knockdown and inhibition led to a similar loss of expression of ß-cell maturity genes in INS-1 cells. Together, these data suggest that STIM1 orchestrates pancreatic ß-cell function and identity through GPER1-mediated estradiol signaling. ARTICLE HIGHLIGHTS: Store-operated Ca2+ entry replenishes endoplasmic reticulum (ER) Ca2+ through reversible gating of plasma membrane Ca2+ channels by the ER Ca2+ sensor, stromal interaction molecule 1 (STIM1). ß-Cell-specific deletion of STIM1 results in a sexually dimorphic phenotype, with ß-cell dysfunction and loss of identity in female but not male mice. Expression of the noncanonical 17-ß estradiol receptor (GPER1) is decreased in islets of female STIM1Δß mice, and modulation of GPER1 levels leads to alterations in expression of ß-cell maturity genes in INS-1 cells.

Pancreatic regional blood flow links the endocrine and exocrine diseases.

An increasing number of studies have demonstrated that disease states of the endocrine or exocrine pancreas aggravate one another, which implies bidirectional blood flow between islets and exocrine cells. However, this is inconsistent with the current model of unidirectional blood flow, which is strictly from islets to exocrine tissues. This conventional model was first proposed in 1932, and it has never to our knowledge been revisited to date. Here, large-scale image capture was used to examine the spatial relationship between islets and blood vessels in the following species: human, monkey, pig, rabbit, ferret, and mouse. While some arterioles passed by or traveled through islets, the majority of islets had no association with them. Islets with direct contact with the arteriole were significantly larger in size and fewer in number than those without contact. Unique to the pancreas, capillaries directly branched out from the arterioles and have been labeled as "small arterioles" in past studies. Overall, the arterioles emerged to feed the pancreas regionally, not specifically targeting individual islets. Vascularizing the pancreas in this way may allow an entire downstream region of islets and acinar cells to be simultaneously exposed to changes in the blood levels of glucose, hormones, and other circulating factors.

Tracking Ca2+ Dynamics in NOD Mouse Islets During Spontaneous Diabetes Development.

The mechanisms accounting for the functional changes of α- and ß-cells over the course of type 1 diabetes (T1D) development are largely unknown. Permitted by our established technology of high spatiotemporal resolution imaging of cytosolic Ca2+ ([Ca2+]c) dynamics on fresh pancreas tissue slices, we tracked the [Ca2+]c dynamic changes, as the assessment of function, in islet α- and ß-cells of female nonobese diabetic (NOD) mice during the development of spontaneous diabetes. We showed that, during the phases of islet inflammation, 8 mmol/L glucose-induced synchronized short [Ca2+]c events in ß-cells were diminished, whereas long [Ca2+]c events were gradually more triggerable at substimulatory 4 and 6 mmol/L glucose. In the islet destruction phase, the synchronized short [Ca2+]c events in a subset of ß-cells resumed at high glucose condition, while the long [Ca2+]c events were significantly elevated already at substimulatory glucose concentrations. In the α-cells, the glucose sensitivity of the [Ca2+]c events persisted throughout the course of T1D development. At the late islet destruction phase, the α-cell [Ca2+]c events exhibited patterns of synchronicity. Our work has uncovered windows of functional recovery in ß-cells and potential α-cells functional synchronization in NOD mice over the course of T1D development. ARTICLE HIGHLIGHTS: In NOD mice ß-cells, 8 mmol/L glucose-induced synchronized short [Ca2+]c events diminish in the early phases of islet inflammation, and long Ca2+ events became more sensitive to substimulatory 4 and 6 mmol/L glucose. In the late islet destruction phase, the synchronized short [Ca2+]c events in a subset of ß-cells resumed at 8 mmol/L glucose, while the long Ca2+ events were significantly elevated at substimulatory glucose concentrations. In the α-cells, the glucose sensitivity of the [Ca2+]c events persisted throughout the course of type 1 diabetes development. α-Cell [Ca2+]c events occasionally synchronize in the islets with severe ß-cell destruction.

Functional characteristics of hub and wave-initiator cells in ß cell networks.

Islets of Langerhans operate as multicellular networks in which several hundred ß cells work in synchrony to produce secretory pulses of insulin, a hormone crucial for controlling metabolic homeostasis. Their collective rhythmic activity is facilitated by gap junctional coupling and affected by their functional heterogeneity, but the details of this robust and coordinated behavior are still not fully understood. Recent advances in multicellular imaging and optogenetic and photopharmacological strategies, as well as in network science, have led to the discovery of specialized ß cell subpopulations that were suggested to critically determine the collective dynamics in the islets. In particular hubs, i.e., ß cells with many functional connections, are believed to significantly enhance communication capacities of the intercellular network and facilitate an efficient spreading of intercellular Ca2+ waves, whereas wave-initiator cells trigger intercellular signals in their cohorts. Here, we determined Ca2+ signaling characteristics of these two ß cell subpopulations and the relationship between them by means of functional multicellular Ca2+ imaging in mouse pancreatic tissue slices in combination with methods of complex network theory. We constructed network layers based on individual Ca2+ waves to identify wave initiators, and functional correlation-based networks to detect hubs. We found that both cell types exhibit a higher-than-average active time under both physiological and supraphysiological glucose concentrations, but also that they differ significantly in many other functional characteristics. Specifically, Ca2+ oscillations in hubs are more regular, and their role appears to be much more stable over time than for initiator cells. Moreover, in contrast to wave initiators, hubs transmit intercellular signals faster than other cells, which implies a stronger intercellular coupling. Our research indicates that hubs and wave-initiator cell subpopulations are both natural features of healthy pancreatic islets, but their functional roles in principle do not overlap and should thus not be considered equal.

High-resolution analysis of the cytosolic Ca2+ events in ß cell collectives in situ.

The release of peptide hormones is predominantly regulated by a transient increase in cytosolic Ca2+ concentration ([Ca2+]c). To trigger exocytosis, Ca2+ ions enter the cytosol from intracellular Ca2+ stores or from the extracellular space. The molecular events of late stages of exocytosis, and their dependence on [Ca2+]c, were extensively described in isolated single cells from various endocrine glands. Notably, less work has been done on endocrine cells in situ to address the heterogeneity of [Ca2+]c events contributing to a collective functional response of a gland. For this, ß cell collectives in a pancreatic islet are particularly well suited as they are the smallest, experimentally manageable functional unit, where [Ca2+]c dynamics can be simultaneously assessed on both cellular and collective level. Here, we measured [Ca2+]c transients across all relevant timescales, from a subsecond to a minute time range, using high-resolution imaging with a low-affinity Ca2+ sensor. We quantified the recordings with a novel computational framework for automatic image segmentation and [Ca2+]c event identification. Our results demonstrate that under physiological conditions the duration of [Ca2+]c events is variable, and segregated into three reproducible modes, subsecond, second, and tens of seconds time range, and are a result of a progressive temporal summation of the shortest events. Using pharmacological tools we show that activation of intracellular Ca2+ receptors is both sufficient and necessary for glucose-dependent [Ca2+]c oscillations in ß cell collectives, and that a subset of [Ca2+]c events could be triggered even in the absence of Ca2+ influx across the plasma membrane. In aggregate, our experimental and analytical platform was able to readily address the involvement of intracellular Ca2+ receptors in shaping the heterogeneity of [Ca2+]c responses in collectives of endocrine cells in situ.NEW & NOTEWORTHY Physiological glucose or ryanodine stimulation of ß cell collectives generates a large number of [Ca2+]c events, which can be rapidly assessed with our newly developed automatic image segmentation and [Ca2+]c event identification pipeline. The event durations segregate into three reproducible modes produced by a progressive temporal summation. Using pharmacological tools, we show that activation of ryanodine intracellular Ca2+ receptors is both sufficient and necessary for glucose-dependent [Ca2+]c oscillations in ß cell collectives.

Exendin-4 affects calcium signalling predominantly during activation and activity of beta cell networks in acute mouse pancreas tissue slices.

Tight control of beta cell stimulus-secretion coupling is crucial for maintaining homeostasis of energy-rich nutrients. While glucose serves as a primary regulator of this process, incretins augment beta cell function, partly by enhancing cytosolic [Ca2+] dynamics. However, the details of how precisely they affect beta cell recruitment during activation, their active time, and functional connectivity during plateau activity, and how they influence beta cell deactivation remain to be described. Performing functional multicellular Ca2+ imaging in acute mouse pancreas tissue slices enabled us to systematically assess the effects of the GLP-1 receptor agonist exendin-4 (Ex-4) simultaneously in many coupled beta cells with high resolution. In otherwise substimulatory glucose, Ex-4 was able to recruit approximately a quarter of beta cells into an active state. Costimulation with Ex-4 and stimulatory glucose shortened the activation delays and accelerated beta cell activation dynamics. More specifically, active time increased faster, and the time required to reach half-maximal activation was effectively halved in the presence of Ex-4. Moreover, the active time and regularity of [Ca2+]IC oscillations increased, especially during the first part of beta cell response. In contrast, subsequent addition of Ex-4 to already active cells did not significantly enhance beta cell activity. Network analyses further confirmed increased connectivity during activation and activity in the presence of Ex-4, with hub cell roles remaining rather stable in both control experiments and experiments with Ex-4. Interestingly, Ex-4 demonstrated a biphasic effect on deactivation, slightly prolonging beta cell activity at physiological concentrations and shortening deactivation delays at supraphysiological concentrations. In sum, costimulation by Ex-4 and glucose increases [Ca2+]IC during beta cell activation and activity, indicating that the effect of incretins may, to an important extent, be explained by enhanced [Ca2+]IC signals. During deactivation, previous incretin stimulation does not critically prolong cellular activity, which corroborates their low risk of hypoglycemia.

Physiological levels of adrenaline fail to stop pancreatic beta cell activity at unphysiologically high glucose levels.

Adrenaline inhibits insulin secretion from pancreatic beta cells to allow an organism to cover immediate energy needs by unlocking internal nutrient reserves. The stimulation of α2-adrenergic receptors on the plasma membrane of beta cells reduces their excitability and insulin secretion mostly through diminished cAMP production and downstream desensitization of late step(s) of exocytotic machinery to cytosolic Ca2+ concentration ([Ca2+]c). In most studies unphysiologically high adrenaline concentrations have been used to evaluate the role of adrenergic stimulation in pancreatic endocrine cells. Here we report the effect of physiological adrenaline levels on [Ca2+]c dynamics in beta cell collectives in mice pancreatic tissue slice preparation. We used confocal microscopy with a high spatial and temporal resolution to evaluate glucose-stimulated [Ca2+]c events and their sensitivity to adrenaline. We investigated glucose concentrations from 8-20 mM to assess the concentration of adrenaline that completely abolishes [Ca2+]c events. We show that 8 mM glucose stimulation of beta cell collectives is readily inhibited by the concentration of adrenaline available under physiological conditions, and that sequent stimulation with 12 mM glucose or forskolin in high nM range overrides this inhibition. Accordingly, 12 mM glucose stimulation required at least an order of magnitude higher adrenaline concentration above the physiological level to inhibit the activity. To conclude, higher glucose concentrations stimulate beta cell activity in a non-linear manner and beyond levels that could be inhibited with physiologically available plasma adrenaline concentration.

pH-Dependence of Glucose-Dependent Activity of Beta Cell Networks in Acute Mouse Pancreatic Tissue Slice.

Extracellular pH has the potential to affect various aspects of the pancreatic beta cell function. To explain this effect, a number of mechanisms was proposed involving both extracellular and intracellular targets and pathways. Here, we focus on reassessing the influence of extracellular pH on glucose-dependent beta cell activation and collective activity in physiological conditions. To this end we employed mouse pancreatic tissue slices to perform high-temporally resolved functional imaging of cytosolic Ca2+ oscillations. We investigated the effect of either physiological H+ excess or depletion on the activation properties as well as on the collective activity of beta cell in an islet. Our results indicate that lowered pH invokes activation of a subset of beta cells in substimulatory glucose concentrations, enhances the average activity of beta cells, and alters the beta cell network properties in an islet. The enhanced average activity of beta cells was determined indirectly utilizing cytosolic Ca2+ imaging, while direct measuring of insulin secretion confirmed that this enhanced activity is accompanied by a higher insulin release. Furthermore, reduced functional connectivity and higher functional segregation at lower pH, both signs of a reduced intercellular communication, do not necessary result in an impaired insulin release.

Dual Mode of Action of Acetylcholine on Cytosolic Calcium Oscillations in Pancreatic Beta and Acinar Cells In Situ.

Cholinergic innervation in the pancreas controls both the release of digestive enzymes to support the intestinal digestion and absorption, as well as insulin release to promote nutrient use in the cells of the body. The effects of muscarinic receptor stimulation are described in detail for endocrine beta cells and exocrine acinar cells separately. Here we describe morphological and functional criteria to separate these two cell types in situ in tissue slices and simultaneously measure their response to ACh stimulation on cytosolic Ca2+ oscillations [Ca2+]c in stimulatory glucose conditions. Our results show that both cell types respond to glucose directly in the concentration range compatible with the glucose transporters they express. The physiological ACh concentration increases the frequency of glucose stimulated [Ca2+]c oscillations in both cell types and synchronizes [Ca2+]c oscillations in acinar cells. The supraphysiological ACh concentration further increases the oscillation frequency on the level of individual beta cells, inhibits the synchronization between these cells, and abolishes oscillatory activity in acinar cells. We discuss possible mechanisms leading to the observed phenomena.

Glucose-dependent activation, activity, and deactivation of beta cell networks in acute mouse pancreas tissue slices.

Many details of glucose-stimulated intracellular calcium changes in ß cells during activation, activity, and deactivation, as well as their concentration-dependence, remain to be analyzed. Classical physiological experiments indicated that in islets, functional differences between individual cells are largely attenuated, but recent findings suggest considerable intercellular heterogeneity, with some cells possibly coordinating the collective responses. To address the above with an emphasis on heterogeneity and describing the relations between classical physiological and functional network properties, we performed functional multicellular calcium imaging in mouse pancreas tissue slices over a wide range of glucose concentrations. During activation, delays to activation of cells and any-cell-to-first-responder delays are shortened, and the sizes of simultaneously responding clusters increased with increasing glucose concentrations. Exactly the opposite characterized deactivation. The frequency of fast calcium oscillations during activity increased with increasing glucose up to 12 mM glucose concentration, beyond which oscillation duration became longer, resulting in a homogenous increase in active time. In terms of functional connectivity, islets progressed from a very segregated network to a single large functional unit with increasing glucose concentration. A comparison between classical physiological and network parameters revealed that the first-responders during activation had longer active times during plateau and the most active cells during the plateau tended to deactivate later. Cells with the most functional connections tended to activate sooner, have longer active times, and deactivate later. Our findings provide a common ground for recent differing views on ß cell heterogeneity and an important baseline for future studies of stimulus-secretion and intercellular coupling.NEW & NOTEWORTHY We assessed concentration-dependence in coupled ß cells, degree of functional heterogeneity, and uncovered possible specialized subpopulations during the different phases of the response to glucose at the level of many individual cells. To this aim, we combined acute mouse pancreas tissue slices with functional multicellular calcium imaging over a wide range from threshold (7 mM) and physiological (8 and 9 mM) to supraphysiological (12 and 16 mM) glucose concentrations, classical physiological, and advanced network analyses.

NMDA receptor inhibition increases, synchronizes, and stabilizes the collective pancreatic beta cell activity: Insights through multilayer network analysis.

NMDA receptors promote repolarization in pancreatic beta cells and thereby reduce glucose-stimulated insulin secretion. Therefore, NMDA receptors are a potential therapeutic target for diabetes. While the mechanism of NMDA receptor inhibition in beta cells is rather well understood at the molecular level, its possible effects on the collective cellular activity have not been addressed to date, even though proper insulin secretion patterns result from well-synchronized beta cell behavior. The latter is enabled by strong intercellular connectivity, which governs propagating calcium waves across the islets and makes the heterogeneous beta cell population work in synchrony. Since a disrupted collective activity is an important and possibly early contributor to impaired insulin secretion and glucose intolerance, it is of utmost importance to understand possible effects of NMDA receptor inhibition on beta cell functional connectivity. To address this issue, we combined confocal functional multicellular calcium imaging in mouse tissue slices with network science approaches. Our results revealed that NMDA receptor inhibition increases, synchronizes, and stabilizes beta cell activity without affecting the velocity or size of calcium waves. To explore intercellular interactions more precisely, we made use of the multilayer network formalism by regarding each calcium wave as an individual network layer, with weighted directed connections portraying the intercellular propagation. NMDA receptor inhibition stabilized both the role of wave initiators and the course of waves. The findings obtained with the experimental antagonist of NMDA receptors, MK-801, were additionally validated with dextrorphan, the active metabolite of the approved drug dextromethorphan, as well as with experiments on NMDA receptor KO mice. In sum, our results provide additional and new evidence for a possible role of NMDA receptor inhibition in treatment of type 2 diabetes and introduce the multilayer network paradigm as a general strategy to examine effects of drugs on connectivity in multicellular systems.

Confocal Laser Scanning Microscopy of Calcium Dynamics in Acute Mouse Pancreatic Tissue Slices.

The acute mouse pancreatic tissue slice is a unique in situ preparation with preserved intercellular communication and tissue architecture that entails significantly fewer preparation-induced changes than isolated islets, acini, ducts, or dispersed cells described in typical in vitro studies. By combining the acute pancreatic tissue slice with live-cell calcium imaging in confocal laser scanning microscopy (CLSM), calcium signals can be studied in a large number of endocrine and exocrine cells simultaneously, with a single-cell or even subcellular resolution. The sensitivity permits the detection of changes and enables the study of intercellular waves and functional connectivity as well as the study of the dependence of physiological responses of cells on their localization within the islet and paracrine relationship with other cells. Finally, from the perspective of animal welfare, recording signals from a large number of cells at a time lowers the number of animals required in experiments, contributing to the 3R-replacement, reduction, and refinement-principle.

ß Cells Operate Collectively to Help Maintain Glucose Homeostasis.

Residing in the islets of Langerhans in the pancreas, ß cells contribute to glucose homeostasis by managing the body's insulin supply. Although it has been acknowledged that healthy ß cells engage in heavy cell-to-cell communication to perform their homeostatic function, the exact role and effects of such communication remain partly understood. We offer a novel, to our knowledge, perspective on the subject in the form of 1) a dynamical network model that faithfully mimics fast calcium oscillations in response to above-threshold glucose stimulation and 2) empirical data analysis that reveals a qualitative shift in the cross-correlation structure of measured signals below and above the threshold glucose concentration. Combined together, these results point to a glucose-induced transition in ß-cell activity thanks to increasing coordination through gap-junctional signaling and paracrine interactions. Our data and the model further suggest how the conservation of entire cell-cell conductance, observed in coupled but not uncoupled ß cells, emerges as a collective phenomenon. An overall implication is that improving the ability to monitor ß-cell signaling should offer means to better understand the pathogenesis of diabetes mellitus.

Corrigendum: Random Matrix Analysis of Ca2+ Signals in ß-Cell Collectives.

[This corrects the article DOI: 10.3389/fphys.2019.01194.].

Random Matrix Analysis of Ca2+ Signals in ß-Cell Collectives.

Even within small organs like pancreatic islets, different endocrine cell types and subtypes form a heterogeneous collective to sense the chemical composition of the extracellular solution and compute an adequate hormonal output. Erroneous cellular processing and hormonal output due to challenged heterogeneity result in various disorders with diabetes mellitus as a flagship metabolic disease. Here we attempt to address the aforementioned functional heterogeneity with comparing pairwise cell-cell cross-correlations obtained from simultaneous measurements of cytosolic calcium responses in hundreds of islet cells in an optical plane to statistical properties of correlations predicted by the random matrix theory (RMT). We find that the bulk of the empirical eigenvalue spectrum is almost completely described by RMT prediction, however, the deviating eigenvalues that exist below and above RMT spectral edges suggest that there are local and extended modes driving the correlations. We also show that empirical nearest neighbor spacing of eigenvalues follows universal RMT properties regardless of glucose stimulation, but that number variance displays clear separation from RMT prediction and can differentiate between empirical spectra obtained under non-stimulated and stimulated conditions. We suggest that RMT approach provides a sensitive tool to assess the functional cell heterogeneity and its effects on the spatio-temporal dynamics of a collective of beta cells in pancreatic islets in physiological resting and stimulatory conditions, beyond the current limitations of molecular and cellular biology.

Heterogeneity and Delayed Activation as Hallmarks of Self-Organization and Criticality in Excitable Tissue.

Self-organized critical dynamics is assumed to be an attractive mode of functioning for several real-life systems and entails an emergent activity in which the extent of observables follows a power-law distribution. The hallmarks of criticality have recently been observed in a plethora of biological systems, including beta cell populations within pancreatic islets of Langerhans. In the present study, we systematically explored the mechanisms that drive the critical and supercritical behavior in networks of coupled beta cells under different circumstances by means of experimental and computational approaches. Experimentally, we employed high-speed functional multicellular calcium imaging of fluorescently labeled acute mouse pancreas tissue slices to record calcium signals in a large number of beta cells simultaneously, and with a high spatiotemporal resolution. Our experimental results revealed that the cellular responses to stimulation with glucose are biphasic and glucose-dependent. Under physiological as well as under supraphysiological levels of stimulation, an initial activation phase was followed by a supercritical plateau phase with a high number of global intercellular calcium waves. However, the activation phase displayed fingerprints of critical behavior under lower stimulation levels, with a progressive recruitment of cells and a power-law distribution of calcium wave sizes. On the other hand, the activation phase provoked by pathophysiologically high glucose concentrations, differed considerably and was more rapid, less continuous, and supercritical. To gain a deeper insight into the experimentally observed complex dynamical patterns, we built up a phenomenological model of coupled excitable cells and explored empirically the model's necessities that ensured a good overlap between computational and experimental results. It turned out that such a good agreement between experimental and computational findings was attained when both heterogeneous and stimulus-dependent time lags, variability in excitability levels, as well as a heterogeneous cell-cell coupling were included into the model. Most importantly, since our phenomenological approach involved only a few parameters, it naturally lends itself not only for determining key mechanisms of self-organized criticality at the tissue level, but also points out various features for comprehensive and realistic modeling of different excitable systems in nature.

Collective Behavior of Social Bots Is Encoded in Their Temporal Twitter Activity.

Computational propaganda deploys social or political bots to try to shape, steer, and manipulate online public discussions and influence decisions. Collective behavior of populations of social bots has not been yet widely studied, although understanding of collective patterns arising from interactions between bots would aid social bot detection. In this study, we show that there are significant differences in collective behavior between population of bots and population of humans as detected from their Twitter activity. Using a large dataset of tweets we have collected during the UK-EU referendum campaign, we separated users into population of bots and population of humans based on the length of sequences of their high-frequency tweeting activity. We show that, while pairwise correlations between users are weak, they co-exist with collective correlated states; however the statistics of correlations and co-spiking probability differ in both populations. Our results demonstrate that populations of social bots and human users in social media exhibit collective properties similar to the ones found in social and biological systems placed near a critical point.

Collective Sensing of ß-Cells Generates the Metabolic Code.

Major part of a pancreatic islet is composed of ß-cells that secrete insulin, a key hormone regulating influx of nutrients into all cells in a vertebrate organism to support nutrition, housekeeping or energy storage. ß-cells constantly communicate with each other using both direct, short-range interactions through gap junctions, and paracrine long-range signaling. However, how these cell interactions shape collective sensing and cell behavior in islets that leads to insulin release is unknown. When stimulated by specific ligands, primarily glucose, ß-cells collectively respond with expression of a series of transient Ca2+ changes on several temporal scales. Here we reanalyze a set of Ca2+ spike trains recorded in acute rodent pancreatic tissue slice under physiological conditions. We found strongly correlated states of co-spiking cells coexisting with mostly weak pairwise correlations widespread across the islet. Furthermore, the collective Ca2+ spiking activity in islet shows on-off intermittency with scaling of spiking amplitudes, and stimulus dependent autoassociative memory features. We use a simple spin glass-like model for the functional network of a ß-cell collective to describe these findings and argue that Ca2+ spike trains produced by collective sensing of ß-cells constitute part of the islet metabolic code that regulates insulin release and limits the islet size.